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I ч 


АВ5ТКАСТ 


The reentry trajectory of a bluff zero-lift reentry body 
was computed using an HP-9845A computer. To obtain the 
external shape of an ablating reentry body as a function of 
altitude, coefficients for polynomial equations were deter- 
mined. At six discrete altitudes, the calculated nose shapes 
were identical to published information for the nose shapes 
of a reentry body with silaceous ablation material. Newtonian 
flow and an exponential atmosphere were assumed throughout. 

The computer program was modified to include ablation due 
to ice crystal erosion, and the trajectory was recalculated 
for a model ice cloud 5,000 feet thick at an altitude of 
30,000 feet. A total of four trajectories was computed for 
various values of mass loss fraction for erosion due to ice. 
The differences in range from the initial trajectory without 


ice erosion were one to fifteen feet. 
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1. INTRODUCTION 


A. BACKGROUND 

The Naval Surface Weapons Center, White Oak, Maryland, 
1s partially responsible for providing improved sensors for 
balloon-carried radiosondes which measure atmospheric proper- 
ties at Various altitudes. Опе item of ongoing research at 
White Oak, and initially the object of this thesis, is an 
improved humidity sensor. 

An object of contention within the strategic missile 
community is whether or not water vapor, rain or ice, i.e., 
atmospheric humidity phenomena, significantly modify the 
reentry trajectory of hypersonic zero-lift reentry bodies. 

In order to determine the influence on the reentry trajectory, 
the ability to measure accurately water vapor, water droplets 
and ice crystals must exist. 

The humidity sensor currently being used on radiosondes 
is a carbon hygristor element which has poor time response 
at high altitudes because of the inherent sensing delays in 
the hygristor element. For temperatures of -40°C, which 
corresponds to altitudes above about 30,000 feet, the time 
lag of the sensing element can be of the order of fifteen 
minutes [Ref. 1]. Therefore, significant errors in humidity 
versus altitude maps are encountered. 

Table 1 is a reprint from Brousaides [Ref. 1] showing the 


time response of the currently used carbon hygristor. Various 
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step changes in relative humidity at several fixed temperatures 
yield the tabulated time responses to achieve 63% or 90% of 


the step input. 


Table 1. Carbon Element Response Lags 


(Reproduced from Brousaides [Ref. 1], p. 16) 


RH Level % Change in RH Indication, sec 
(Initial-Final) Temp 3 90% 








О) 
егі 





30-50 0.1 0.5 | 
51-85 0.3 1.4 

85-51 0.4 2.2 

51-30 0.3 1.6 | 
81-21 2.6 14 

50-80 2,4 21 

79-51 2.6 48 

51231 3. 8 41 

31-50 15 61 

50-72 11 104 

72-50 9 169 

50-31 33 163 

31-50 86 363 

50-66 82 212 

55-50 88 578 

50-30 260 727 


Table 1 shows that, as temperature decreases, that is, as 
altitude increases, the time response of the element becomes 
so slow as to make the readings virtually meaningless for 
any reasonable rate of upward balloon motion. Іп the case 
where ice is to be detected, the time lag of the sensing 
Element is particularly critical. 

The object of current research on the hygristor element 
is to determine a heat transfer model of the sensing element 
so that the delay times can be calculated as a function of 


temperature. By relating delay time to temperature and 
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temperature to altitude, a more accurate map of humidity versus 
altitude can be obtained. If an accurate heat transfer model 
of the hygristor element can be obtained, data processing can 
be used to yield accurate measurements. Alternatively, an 
improved sensor can be designed once the exact cause of the 
delay can be pinpointed. 

Even if humidity can be accurately measured, determination 
of water droplet or ice crystal concentration and size distri- 
bution within a cloud is difficult. In other words, even if 
humidity can be sensed accurately, can this information be 
related directly to rain or ice? Perhaps a different type 
of sensor or several sensors might be needed to determine 
humidity, cloud composition and concentration. 

In addition, in order to design a useful sensor, criteria 
for measurement accuracy are required. The criteria for 
sensor accuracy are determined by the sensitivity of the 
reentry trajectory to changes in water vapor, rain and ice 
concentrations. 

Criteria are not available for determining the change 
in a reentry trajectory as a function of atmospheric water 
content. Therefore, it became necessary to start at the 
beginning in order to learn enough about the reentry process 
to make estimates of how much the reentry trajectory could 
differ with changes in humidity-related atmospheric conditions. 
This thesis estimates the change in range due to clouds with 


wee crystals. 
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Consequently, basic information concerning hypersonic 
flow about bluff bodies was studied. Equations of motion 
for reentry bodies were formulated in relation to the current 
problem. The many factors relating to reentry, such as 
ablation, convective and radiative heating and the modifi- 
cations to ablation because of chemical interaction with 


water were considered. 


гэ SCOPE 

Three separate atmospheric water phenomena affect 
reentry ablation. The three phenomena are water vapor, rain 
and ice crystals. ОЁ the three phenomena, reentry through 
ice clouds is the most destructive. Consequently, this 


thesis studies the effects of ice clouds. 
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JI. THEORY 


A. REENTRY DYNAMICS 
ШО Trajectory 

The trajectory of interest is the ballistic reentry 
of a zero-lift vehicle typical of reentry bodies used as 
warheads. A ballistic missile trajectory is typified by 
near-orbital reentry velocities and a steep initial flight 
path angle, i.e., greater than five degrees. 

lI! trajectories Of reentry bodies pass through 
three operational phases. The three operational phases are 
defined as follows: 

a. Keplerian phase: The segment of the trajectory 
where gas dynamic forces are negligible (p = 0). 

b. Intermediate phase: The segment of the 
trajectory where gas dynamic forces are of comparable 
magnitude to other terms (gravity) in the equations of 
motion. 

c. Gas-dynamic phase: Тһе segment of the trajectory 
where aerodynamic forces dominate the equations of motion 
for reentry. 

Of the three phases described above, the gas dynamic 
phase is the one in which environmental factors, such as 
clouds or humidity, influence the trajectory. Therefore, 
the gas dynamic phase is the only portion of the reentry 


which will be considered. 
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The analysis of the flight of the reentry body is a 
particular application of the general theory of the dynamics 
of rigid bodies in three dimensions. The motion of the center 
of mass can be separated from the motion of the body about the 
center of mass. In this discussion, since a zero-lift body 
is hypothesized, the assumption is made that the center-line 
of the body always is parallel to the velocity vector. 

In addition, the assumption will be made that the 
reentry trajectory is planar rather than three-dimensional. 
Motion in the third dimension is, of course, possible and 
indeed certain for actual reentry bodies due to wind, 
precession, Coriolis force and so forth. However, if the 
"particle" can be described as one capable of generating drag 
only, then motion in only two dimensions is assured.  Conse- 
quently, the reentry body appears as a point particle of 
varying mass subject to the various forces applied by the 
entry environment. 

Atmospheric properties, particularly the densities 
existing at higher altitudes, are subject to considerable 
uncertainty. In order to provide a model of the density 
profile of the atmosphere, two assumptions are made: 


l. The earth and its atmosphere are spherically 
symmetrical. 


2. The atmospheric density varies exponentially 
with altitude. 


The first assumption is reasonable for this analysis 
since only a small portion of the atmosphere is traversed 


for a ballistic entry. The second assumption is based 
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on the simple kinetic theory of an isothermal gas in a 
uniform gravitational field. These assumptions lead to the 


well-known exponential model for the atmosphere: 


p = p e m (1) 
where: o - free stream atmospheric density (slug/ft°) 
2 assumed Mean sed level atmospheric density 
B = atmospheric density decay parameter (ft +) 
y = altitude of vehicle above the surface of 


the earth (ft) 
B is equal to the mean molecular weight of the earth's 
atmosphere times the local acceleration due to gravity 
divided by the product of the mean atmospheric temperature 
and the universal gas constant. Figure 1 1S a reprint from 
Loh [Ref. 2] and compares the exponential approximation 
with the ARDC (1956) model of the atmosphere. 
2. Equations of Motion 
The geometry of an IRBM or ICBM trajectory is shown 
in Figure 2 [Ref. 3]. Assuming no friction occurs during 


the flight, the conservation of energy gives: 


2 


1/2 mV? - mer = 1/2 nV, 


- MgR = constant (2) 


Equation (2) provides a unique relationship between 


velocity and altitude provided some initial condition 
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Figure 1. Comparison of Exponential Approximation 
with ARDC Model of Earth's Atmosphere (1956) 
(Reproduced from Loh [Ref. 2], p. 13) 





Figure 2. Diagram of ICBM Flight Path 
(Reproduced from Fuhs [Ref. 3]) 
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or boundary condition can be used to find the constant term. 


~g 


Equation (2) also can be written, if r “ В, ав: 


Rm 


ат? 
SS = 


m 
7 ас 


(г 9992] - g = constant (5) 


where the constant in both equations (2) апа (3) is the total 
reentry energy, i.e., kinetic plus potential energy, of the 
body. 

The conservation of angular momentum of a constant 


mass also applies and can be written as: 


mr w = тг“ = = mRV, (0) = mRV ,cosA (4) 


Boundary conditions for equations (3) and (4) ате 9(0) - 99 
and r(0) = R. Time can be eliminated from the equations 


using the relation from above. 


dr атсаз dr no 


19 z- = 
T. 


Substitution of equations (4) and (5) into equation 


(3) gives: 








ae 2 
- 25 2 









21 





Let 


RKV COSA 


И ф _ gR 
ON — V,cosal (7) 
Note that 
-RV,cosA 
dZ = Шинээ dr (8) 
T 


Substituting equations (7) and (8) into equation (6) yields 





(9) 
(10) 
E: 
о к i 11 
ф Б IGE (11) 


then equation (10) becomes 
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тг. e) 


ф 


кіре 


To evaluate the constant C, symmetry and the boundary 


conditions r = R, Ө = АР апа 2 = 2, сап БЕ ИСО : 
А R 
Z со 703 | 
ne ф 
т s. — 3 
gR _ 28R 
| . = V : 
ф ф 
Substitution of the boundary conditions into equation (12) 
pads: 
ES 
5 _ 2 oL (E4) 
COS (9. С) cos (8, (БЛ ф 
This implies that C = m and therefore: 
_ m 
cos (0 + T) = - cos = $ (2125 ) 


At this point, refer to equations (2) and (3), and 
note that for some value of initial velocity the total energy, 
i.e., the constant term, will be zero. This value of initial 


velocity corresponds to the initial velocity at which the 
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particle kinetic energy is equal to potential energy. Since 
the constant term is zero, as the particle moves, the constant 
term will remain zero, and kinetic energy will just balance 
the potential energy for all values of r, even as r approaches 
infinity. Thus, the value of initial velocity is the escape 


locity. Hence 


V = V “ = 208 (16) 


Substitute equations (7) and (11) into equation (15) to 


get 
gR се 
V„cosA T 
ОО Вэ  ----- (17) 
2 
2 gR 
% 3 2-2 - 2gR 

Define and substitute the relations V - VAR and r' = r/R 


and equation (16) into equation (17) to arrive finally at 


2 r (18) 


cosg = 5 
(1 + 4V cos A(V  - 1)) 


Equation (18) is used to calculate a family of curves 


showing range versus launch angle. Іп Figure 3, range is 
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Bund as a function of launch angle be setting r' = 1,0 in 
equation (18) to find m for various values of V between 


0.5 and 0.7 and realizing that 


Бапве - 280. (19) 


where R = 2.09 x 10% E Tor the radius of the earth. 


“2” 256074 ft/sec 


т «0.55 
Fg = 23310 Зэс 


Y = 0.50 
vg * 22005 <t/see 


Y = 0.55 
v3 * 20172 2t/see 


y e 05 
77” 18335 2*/sa- 


RAUCH, THOUSAINS OF NAUTICAL MILES 





9 20 20 70 40 50 50 79. 89 50 
LATNCS AXGLZ 


р 


Figure 3. Range of a Ballistic Missile as a Function of 
Launch Angle for Various Values of V = V al Mes 
(Reproduced from Fuhs [Ref. 3]) 
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A value of e can be deduced from Figure 3 by assuming 
that the missile will be launched with the proper velocity 
and launch angle to achieve the optimum range, i.e., the peaks 
Of the curves. 

For example, Table 2 shows the various characteristics 
of typical ballistic missiles in the inventory today. The C-4 
Trident has a range of between 4000 and 4800 nm. Referring to 


Figure 4, an optimum range of 4100 nautical miles requires a 


launch angle of about 27? and 1 of about 22006 ft/sec. 


Table 2. Launch Vehicles and Parameters 
(Reproduced from Fuhs [Ref. 3]) 



























Misstle da 


ЖАЛДАП 


weight Payload 





Pange Warhead 106 


run Yield 


Retired 























А-1 Folaris 30,000 Мау 60 Oct 65 2 
A-2 Polaris 54"31' EE Ea 
A-3 Polaris 5421! 2 
C-3 Poseldon 7]2"34' 1800 

C-4 Tricent 72" 31! 70,900 1500 3 
0-4 Trident 82"45,8! 

0-5 Trident M 126,000 


MX 227 


Atlas 


00" due to 
guidance) 
30' 120) 
A шиг. | 
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Returning to equation (2), the constant for the C-4 


Trident can be found by substituting V, and R. 


Z 


) constant _ x (20) 


- gR = л 


I 
7 V 


8 2 


K = 5 (22006)° - (-32.2)(2.09 x 107) = 9.15 x 10% £t”/sec 


Since K now is known, equation (2) can be used to find the 
ENSerty at any point in the trajectory given r by rewriting 


equation (2) as 


у = [2(K + gr)]1/? (21) 


If the portion of the trajectory near reentry is 
examined (that is, where r is not much different from R), 
the earth can be approximated as being flat since y/R << 1 
where y is altitude. For this case, the simple two-dimensional 
equations of motion in Cartesian coordinates apply, that is 
2 mx 222) 


апа 


ХЕ = my (23) 
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where x is distance measured horizontally, and y is distance 
measured Е слу. 

Since in the model the only force in the x-direction 
1s drag, and the forces in the y-direction are drag and 


gravity, equations (22) and (23) can be written as 


De 

> m 222) 
р 

Мы q^ zo (25) 


where pe and D, Lepre Seneste drag in the x and y directions, 
respectively. In finite difference form, equations (24) and 


(25) are 


2X ; 1 X; Dy 
т (26) 
A 
А шоуг р 
+ 
2 l= -g - Z (27) 
t 


The values of X1 and Yi and V(0) are obtained from initial 
conditions at the point of interest at the beginning of 
reentry using equations (18), (19), and (21). The values of 


X, and Y, are obtained from 





aa, Е U 5t (28) 


< 
1 


е Vent (29) 


where 


(30) 


ШЕ ПҮРЕК5ОМІС КЕЕКТЕҮ 

The character of the gas dynamic phase of flight at 
near-orbital velocities is far removed from that of flight at 
sonic or subsonic velocities. The vehicle has tremendous 
energy and angular momentum, most of which must be released 
to the atmosphere. The two processes available to transfer 
the energy are manifested as heating of the vehicle and 
deceleration. 

No precise definition may be given for the velocity 
at which supersonic flow becomes hypersonic flow because 
the effects vary with the vehicle geometry. In addition, the 
onset of these effects occurs gradually rather than abruptly. 
Knowledge of Mach number at which transition from supersonic 
to hypersonic flow occurs is important. If the flow is hyper- 
sonic, certain simplifying assumptions can be made concerning 
the shock shape and position relative to the body, which 
relate directly to drag and pressure coefficient. These 


assumptions are discussed below, 
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During hypersonic flight, there will be regions of the 
flow fields about a typically blunt-nosed body where the 
flow decelerates and in which subsonic, transonic and 
supersonic flows all exist simultaneously. For the purposes 
of this discussion, a Mach number greater than or equal to 
four will be considered hypersonic. 

Any bluntness of the nose will manifest itself as a 
detached shock. A pointed nose will melt or ablate due to 
the intense heat generated at the point and therefore will 
become blunt. Consequently, the blunted nose generates the 
detached shock. The most striking effect of hypersonic flight 
is the formation of a shock ahead of virtually all bodies. 
The behavior of the shock for the flow past a blunt body of 
revolution as the Mach number increases may be illustrated 
as in Figure 4, taken from Cox and Crabtree [Ref. 4], for 
flow of a perfect gas about a sphere. 

As can be seen, the shock stands off from the body at all 
Mach numbers, and the shock-detachment distance decreases 
as the Mach number increases. The limiting shock (M > œ) lies 
close to the body and is very nearly concentric over most of 
the surface of the sphere. Just behind the shock wave where 
the shock is nearly perpendicular, the flow is brought almost 
to rest and is subsonic. As the flow passes around the body, 
the flow accelerates, passes through sonic speed and finally 
accelerates to supersonic velocities as is illustrated in 


Lure 5, again taken from Cox and Crabtree [Ref. 4]. 
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Figure 4. Shock Shapes for Flow Past a Sphere at Different 
Mach Numbers (Y = 1.4). (Reproduced from Cox 
and Crabtree [Ref. 4], p. 9) 


At the high temperatures experienced behind the shock in 
hypersonic flow, the assumption of a perfect gas cannot be 
made as the energies of the gas associated with temperature 
are comparable to the energies of various molecular and atomic 
processes such as ionization, dissociation and vibration. 
Therefore, characteristics of the gas such as the commonly 
accepted value of the ratio of specific heats for air, l.e., 

Y = 1.4, must be viewed with skepticism and taken into account 
Men calculating the flow field. 

However, the Mach number can be associated with the ratio 

of the kinetic energy and the random thermal energy by the 


relation 


S 
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subsonic and Supersonic Regions for Flow Past a 
Blunt Body. (Reproduced from Cox and Crabtree 
| Ке ТШ, р 65) 


large Mach numbers, regardless of the value of y, the 
energy is much larger than the thermal energy. This 
that the random motion associated with collisions 


atoms and molecules is of relatively little importance 


compared to the kinetic energy of the particles and therefore 


the gas may be thought of as a stream of non-interacting 


particle 


5. The stream of non-interacting particles was first 
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conceived by Newton for a "rare" gas composed of particles with 
large distances between them. Newton showed that if the 
particies do not interact with one another and at large 
velocities, then the resistance of a body into which the 
particles collide will vary as the square of the velocity. 

In addition, he calculated, by assuming that there was no 
transfer of tangential momentum and that the normal momentum 
of the particles is destroyed at the surface, that the force 
on a plate at incidence to the flow varies as the square of 
the angle of incidence. More recently these assumptions have 
been used to show that the pressure on the plate is given by 


AP = DV siny (52) 


от 


28110 (35) 


©? 
ii 


Also, for flow over wedges at high Mach number, the 
assumptions were used to show that the relation between the 


shock angle (u) and the body angle (v) is 


(34) 





Therefore as y > 1, the body and the shock coincide and 
the density ratio across the shock becomes infinite. Thus 


there is a layer of high density gas lying near the surface. 
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This is the "Newtonian Shock Layer" or entropy layer since 
very large changes in entropy also occur. 

Lees [Ref. 5] adjusted the Newtonian Formula [equation (33)] 
empirically so that the value at the stagnation point on the 
nose is exact. Lees' result, which is called the "Modified 


Newtonian Formula," is given by 


С = C sing ОС; 


where С, is the stagnation point pressure coefficient. Тһе 
О 


value of Cy is given by the Rayleigh formula developed from 
О 
the integral Bernoulli equations and the shock equations and 


1s 


1 Y 


p 2 2 2 
О ҮМ 2M y + 1 


As in the Newtonian formula, as M > œ and y > 1, the value 








mt C approaches two. 
Po 


There are various other modifications to the Newtonian 
Formula including the Busemann Formula, the Modified 
Newtonian-Busemann Formula and others but, as is shown in 
Figures 6 and 7, a comparison of the surface pressure 
coefficients for the flows over a sphere and a cylinder at 


M > «e and y = 1.4 confirms that the Modified Newtonian 
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Figure 6. Surface Pressure Distributions for Flow Past a 
Sphere (M = =, y = 1.4). 
(Reproduced from Cox and Crabtree [Ref. 4], 
p. 75). 
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Figure 7. Surface Pressure Distributions for Flow Past a 
Circular Cylinder (M = », y = 1.4). 
(Reproduced from Cox and Crabtree [Ref. 4], 

ВЭ 1241 


36 





Formula is in best agreement with the actual numerical 
calculations. 
1.  Aerodynamic Heating and Ablation 

During hypersonic reentry, excessive local heating 
15 а serious problem. For the ballistic vehicle, where the 
lift-to-drag ratio is zero, very high heat fluxes are 
encountered but for a relatively short duration, i.e., five 
to 100 seconds. In practice, the heat and mass transfer 
phenomena associated with this reentry mode are nonlinear 
functions of time due to the nonlinear environmental parameters, 
chemical reactivity, mechanical and acoustical vibration, 
gas temperature, mode of energy transport, total heat load, 
peak heating rate, shape of the heat pulse, particle impinge- 
meat, Material properties, etc. Of primary concern is the 
means of accommodating the energy converted to heat. 

The thermal protection of the reentry body can be 
accomplished in three ways: (a) absorption, (b) radiation, 
and (c) reduction of aerodynamic heating prior to reaching 
the vehicle. The added requirement to minimize body weight 
in order to achieve greater missile range has often led to 
the choice of materials designed to change phase when heated 
and to vaporize at the surface exposed to the hypersonic 
environment. The energy absorption achieved through such 
processes has proved to be highly effective. 

The energy transfer mechanisms between the boundary 
layer gas and the ablating surface are shown in Figure 8, 


Beprinted from Hurwicz, Kratsch, and Rogan [Ref. 6]. The 
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Figure 8. Heat Flux Balance at the Surface. 
(Reproduced from Hurwicz, Kratsch, and Rogan [Ref. 6], p. 9) 


energy balance at the surface established the boundary con- 
dition for solution of the boundary layer equations. The net 
heat flux through the surface to the interior of the ablator 
(q,) equals the summation of the aerodynamic heat flux to 

the surface from the boundary layer (ae the energy trans- 
fer away from the surface by normal gas mass flow А у the 
absorption of incident radiation C> the reradiation of 


4 


energy from the surface (te, oT, ), the energy transfer by 


mechanical erosion (ablation) (Ah), the energy transfer to 


the surface by condensed material (m. h.), and the energy 
S ~s 


transfer to the surface by internal gas products m, h, Dl 
5 95 


Note that the aerodynamic heat flux am is the combination 
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of convection to the surface as well as chemical energy 
Master due to surface chemical reactions from mass transfer. 


The preceding factors can be expressed as an energy equation 


et, mh tm h. + "gg. 122 mh. (97) 


Equation (37) is applicable for both laminar and 
turbulent flow and couples the energy balance to the boundary 
layer solution which is required to obtain species concentra- 
tions and the heat and mass transfer analysis. 

Equation (37), although general, does assume thermo- 
chemical equilibrium at the surface and in the boundary 
layer. Such equilibrium may not exist in reality since the 
mass ablation rate probably is reaction rate controlled by 
the ablation material, its char products, and the boundary 
layer species. In this case, gaseous ablation products can 
be determined in principle through the simultaneous solution 
of: (a) all surface thermochemical equilibrium reactions 
expressed in partial pressures and equilibrium constants, 

(b) equations of conservation of elemental species at the 
surface, and (c) Dalton's law of partial pressures. 

Even though conceptually the surface thermochemical 
equilibrium reactions can be expressed explicitly, the actual 
products of the reactions are difficult to determine because 


of the elevated temperatures and the number of interacting 





elements. If, for example, carbon is used as an ablator, 
me products normally considered are C, C» Ez» Сі, Cc, СО, 


20 EN, О М,. As can be seen from Figure 9, the concen- 


pr ==? 

trations of the gases change as temperature (or velocity) 
changes. In addition, if water vapor is included, as is of 
interest in this discussion, another element, hydrogen, must 


be included as a possible reactant. Thus, entire organic 


chains of molecules may be formed. 


0 E 
“> 
10° ae 
<a 
6 | 
nt m & | | 
10° A : | 


0 10,000 20,000 30,000 
VELOCITY FT/SEC 


CONCENTRATIONS 
i 
ё 


Figure 9. Mole Fraction Concentration behind Normal Shock 
(Reproduced from Cox and Crabtree [Ref. 4], p. 10) 


In view of the complexity of ablation with non- 
equilibrium chemical compositions, information concerning 
nose shape was sought in lieu of direct calculation.  Informa- 
tion was found relating nose shape to specific altitudes. The 


nose tip chosen was composed of a silaceous material. 
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Clon occurred during reentry through the atmosphere 
including mechanical erosion. Further, turbulent flow was 
assumed to exist from the start of reentry and continue 
throughout the flight. Figure 10 shows the family of nose 
shapes selected for this analysis. The material loss due to 
ablation is excessive compared to modern carbon nose tips; 
however, the effects of nose shape change on the trajectory 
were to be examined. Therefore, the larger the ablation 
rate, the more the trajectory could be expected to change 


from the non-ablating entry. 
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Figure 10. Nose Shape of Reentry Body at Various Altitudes. 
(Reproduced from Hurwitz, Kratsch, and Rogan 
езара i 


Recent experiments have shown that interesting after- 


body ablation patterns are formed during reentry. 
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Gross-hatching or other similar spiral ablation patterns 
occur, and related roll torques are induced producing erratic 
perturbations in roll behavior. Such phenomena were not 
considered in this analysis. 
2. Drag 
The drag of a reentry body experiencing Newtonian 


flow can be expressed as: 


Drag = q C sin) dA (58) 


پس 


where q represents dynamic pressure, C, represents pressure 


p 
coefficient, ф 15 the local slope of the body, and A is an 
element of area. Given the nose shapes in Figure 10 and 
assuming that the shape of the nose can be determined for 
every altitude, the drag at any altitude can be computed 
using equations (35) and (36). For an axisymmetric body 
where т - f(x), the body сап be divided into infinitesimal 
Segments of X; and the drag can be found on each elemental 
circular element. Figure 11 shows a typical element in the 
X-r plane. Figure 12 shows the same element as seen along 


the longitudinal axis and where R is the base radius of the 


reentry body. Therefore, the incremental area is 


_ EE UD dr 
dA 2Trds "Sind (39) 
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Figure 11. Cross-Section of Elemental Area along the Nose. 


Figure 12. Front View of Elemental Area along the Nose. 
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Hence, equation (38) can be written for the elemental drag 


and becomes 


Y tr 
d(Drag) = ЗЕ лга? (40) 
Т 


Tbe total drag on the body is then the sum of the 


elemental drags over the length. 
2 


Since the Newtonian flow model applies over the 
length of the nose, the 'afterbody" can be included in the 
nose shape equations if properly written; hence equations (38) 
to (40) will define the drag in all regions along the body. 

For projectiles which fly at low supersonic velocities, 
the base drag 1s very important when considering the total 
drag on the body. However, at hypersonic velocities, the 
pressure increase across the shock at the nose is so large, 
of the order of one hundred atmospheres, that the low pressure 
at the 2 the reentry body is comparatively insignificant. 


In this analysis, the base pressure drag was neglected. 
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See HUMIDITY EFFECTS 

As noted previously, the effects of water vapor, water 
ets and ice crystals іп the atmosphere lead to extremely 
complex modifications to all the equations previously cited. 
Considering water vapor only, density of the air is changed. 
As a result, the modified density influences reentry body 
velocity, ablation rate, etc. Also, as mentioned before, more 
complex chemical equations are required to determine ablation 
products and, hence, heat loss and absorption. The total 
effect on drag, for instance, would be the sum of the effects 


on each parameter affecting drag; for example: 


Das — A denssty . A E qo Mna (42) 
Drag density Ср q area 2 


If rain, i.e., suspended water droplets, also is included 
in the analysis, mechanical erosion must be considered. 
Additionally, when the droplets cross the shock wave in front 
of the vehicle, scattering and annihilation will occur; and 
only some fraction of the droplet mass actually will contact 
Me Surface of the body. The interaction of water vapor and 
raln with the reentry body is not considered in this thesis. 

The interaction of ice in high altitude clouds would seem 
to be a worse case of humidity effect interaction with a reentry 
vehicle since the assumption can be made that most of the 


typical ice particles will survive after passing through the 
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shock front. Hence, mechanical erosion should be increased 
greatly. Experiments have shown that the mass ablated can 
be many times the impact mass of the typical ice particle. 

Equations also are simplified since the time in the 
boundary layer is so short, of the order of microseconds, 
that the ice will not undergo a change in phase before 
colliding with the body. Therefore, changes in ice crystal 
EESDerties and chemical interactions were neglected. 

Since information was not available concerning the 
actual change in nose shape due to collisions with ice 
particles, the assumption was made that ablation shape was 
not changed but merely was increased in rate of mass loss. 
That is, the identical equations describing the nose shape 
were used in both the baseline trajectory model and the ice 
cloud trajectory model. However, as mass loss increased 
due to impact, the nose shape was changed according to those 
equations until the proper mass loss was achieved. This 
assumption is no doubt incorrect in that it would be expected 
that the mass loss would be much greater near the stagnation 


point and much less along the perimeter of the body. 
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III. SIMULATION 


Pee BASELINE TRAJECTORY 

In order to evaluate the consequences of reentry through 
ісе clouds, a reference trajectory was needed to compare such 
parameters as range, impact velocity, changes in drag coeffi- 
cient, etc. The equations of motion were presented previously, 
EwrEusnitial conditions and boundary conditions are required 
to solve for a specific trajectory. 

1. Selection of a Typical Reentry Vehicle 

The model selected as the reentry vehicle conforms 
to the nose shapes shown in Figure 10. In addition, realistic 
values of ballistic coefficient and drag coefficient were 
required. Modern ballistic reentry vehicles have ballistic 
coefficients (weight divided by the quantity drag coefficient 
times base area) in the range of 1500-2000 pounds per square 
moot and drag coefficients slightly greater than 0.1. 
The weight of the reentry vehicle was not known, 

but a specific gravity of 1.4 was selected as typical. 
During the course of the initial iterations of the computer 
model, the body length was changed to add enough mass so that 
@ac ballistic coefficient was near the desired value of 2000. 
The body shape was achieved by extending the after-section of 
the initial nose shape shown in Figure 10. Figure 13 shows 


the final model shape prior to reentry. 
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E initial Conditions 

Nerterteina reo rigure 4, the center trajectory was 
chosen as a representative trajectory of a typical Navy 
reentry vehicle which might be launched by Polaris or Trident. 
This curve corresponds to a maximum range of about 4100 
nautical miles if launched at an angle of 27° with an initial 
velocity of 22006 feet per second. Substituting the values 
for launch angle and initial velocity, equation (20) yields 


E £t^/sec* as shown previously. 


9.15 x 10 
Since deceleration and ablation have not yet begun 
at 250,000 feet, this altitude was chosen as the initial 
altitude to begin reentry. For an altitude of 250,000 feet, 
the value of r, the distance to the center of the earth, is 


20115 х 10” feet; the value for r was substituted into 


equation (21) to yield 


Vi = [2(K + 2315/2 огото х 10% - 32.202.115 x 10971? 


21652 tSc (43) 


At an altitude of 250,000 feet atmospheric density is too low 
to cause drag. Also, at that altitude r is not much different 
from R, the radius of the earth. Therefore, the coordinate 
system can be shifted from r,9 to x,y. Since x forces are 


zero, the x component of velocity will remain constant and 
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ШЕШ! De equal to the launch velocity times the cosine of the 


launch angle 


Wie = V ¿COSA >. 2006 EOS 27° = 193607 £t/sec (44) 


At 250,000 feet, the x-component of velocity must be 
the same value; and since the value of V is known, the 


reentry angle À can be determined from 





A = cos”? c = cos Í 12907 = 25.04° (2125) 
ф 


Checking these values with the curve selected in 


Figure 10 and using equation (18), 


E 2 o 
Е. - О) (сос (272771 = 0.827 


[1 + 4(.06)% cos*(27%)(.6% - 1)] 1/2 


Therefore, the initial angular distance from apogee is 


E 34.21%. Thus total range using equation (19) is 


7 
то 0220:1221)/ї- 


Range = 2R9 
5 5 5000 


which corresponds very well to the expected value of 4100 nm. 
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The range the reentry vehicle is expected to fly сап 
be found from the reentry location to impact. То find the 
Ie for a no-drag trajectory recall that the arc length 


traveled is 
0! 20 - 9 (46) 


Again, using equation (18) and substituting values for 


I and r' at the entry point 


DEEP 2 Š 
OBB n --022::0127/(1101) мэ. 


[1 * 4(.6)^cos^(25.04?) (.6^ - 1)] 


and 9 = cos *(.843) = 32.54%. 


Miewangular arc of flight using equation (46) is 
00-9: 51219 - 52,549 = 1.67° 
and using equation (19) 


7 
" NP o POSEE 077. 
Range R9 00 a Е OPES sala 
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Since the sequence of equations above assumed no drag, the 
range should be the maximum expected range of the model 


reentry vehicle. 


Table 3 


Summary of Initial Conditions at 
an “Аш итде of 250.000 Feet 


Initial Velocity 21,632 ft/sec! 
Initial Entry Angle 25° | 
Initial x velocity 195607 ft/sec 
Initial y velocity 9,138 ft/sec 
Maximum range OS nm | 


3. Nose Shape 
The nose shape is the primary factor in determining 

total vehicle drag. Since the nose shape was given only for 
discrete values of altitude as in Figure 10, the nose shape 
WAS constructed artificially as a continuous function of 
altitude. This was done by using the НР-9845А computer along 
with the HP-9874A Digitizer. Each of the six nose shapes was 
divided into two segments; one segment was from the nose tip 
to 2.4 inches aft of the nose tip, and one segment was from 


Mee inches aft of the nose tip to about 20 inches aft of the 
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nose tip. Another function was required to describe the 
recession of the nose tip. 

Each nose shape was traced with the digitizer 
inputting data points into the computer. A curve fitting 
routine was used to construct curves of nose radius as a 
function of length along the centerline of the reentry body 
using the data from the digitizer. The best fit curve for 
the segments from the nose tip to 2.4 inches aft of the nose 
tip was a fifth-order polynomial with constant coefficients. 
The segment from 2.4 inchés aft of the nose tip was described 
by a straight line. Appendix A contains graphs of each of 
the six nose shapes (two segments each) comparing the data 
input with the plot of the equations which describe the data. 
As can be seen, the graphs very closely approximate the 
actual nose shapes. 


The nose shapes thus were defined mathematically as 


5 
1 
r = ¢f А.х 
i=0 1 (47) 
mero < x < 2.4 inches and as 
1 i 
r= 5 Bix (48) 
1=0 
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for x > 2.4 inches where r is the local radius of the 

reentry body and x is the distance along the centerline. 
Using the coefficients calculated in equations (47) 

and (48), plots of each coefficient versus altitude were 

E ructed. For example, all the fifth-order coefficients 

of equation (47) were plotted at the altitudes applicable 

to the particular nose shapes. The coefficients in equation 


(47) thus were defined by 
6 ! 
шаг (49) 


where h is altitude in feet. The coefficients in equation 


(48) were defined by 
6 Р 
B. = < OD..h? (50) 


Sion of the coefficient curves into two portions was 
necessary; equations (49) and (50), sixth-order polynomials, 


applied below 90,000 feet. Equations (51) and (52) 
1 
A. = у E..h? (51) 
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73) (52) 


were used above 90,000 feet. 

Equations (47) through (52) can be combined to 
N the nose shape, the radius, as a function of altitude 
in the various regions of applicability. Table 4 shows the 
combined equations (53) through (56) and the regions of 


Eusrcability. 


Table 4 


summary of Nose Shape Equations 


ACES 


x > 2.4 (inches) 


Equal to Б 6 


| 
| 
. | 1 6 EV 

= E 2 C,;h’x' (53) тег D, hx (54) 
feet Ln | KO Jag 

| 
90,000 E AS IO y Ph xe (56) 
feet i=0 j=0 12 150 1501 


Appendix B compares graphs of the polynomials 
@escribing the А, coefficients versus altitude with corre- 
Spending graphs of the digitizer data input used. Appendix C 
SemtainSs Similar graphs of the nose tip recession versus 
altitude. All graphs show that the equations adequately 


represent the described data. 
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The result of the above analysis was to allow the 
insertion of any altitude inco cquations (53) through (56), 
yielding the corresponding nose radius as a function of the 
Centerline distance from the nose tip. 

4. Flow Chart 

Using the initial entry conditions and nose shape 
versus altitude information previously described along with 
the equations of motion, a program to calculate trajectory 
was possible. Figure 14 describes the process. 

5. Validity 

The program initially was modified to print the nose 
shapes every second along the trajectory so that a comparison 
between the desired nose shapes described in Figure 10 and 
the nose shapes described by the equations derived in 
Section IIIAS could be made. Figure 15 shows a family of 
eurves plotted from the program output for nose shape at 
altitudes closely corresponding to the altitudes for the family 
of curves in Figure 10. Although the nose shapes in Figure 15 
are somewhat more square in appearance than those in Figure 10, 
the proper general shapes were clearly achieved with the nose 
shape as a function of altitude equations. 

A goal of the model was to achieve a drag coefficient 
of about 0.1 and a ballistic coefficit of about 2000. In the 
simulation, the drag coefficient ranged from 0.100 to 0.127 
during the flight with the initial value of 0.105. The 
ballistic coefficient ranged from 1392 to 1878 during the 


flight with the initial value of 1798. The range should be 
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ios than the vacuum trajectory range of 101.5 nm; the range 
with drag was 84.103 nm. 

TOT ballistic entries at escape velocity and initial 
angles of inclination greater than about 5°, the variation of 
yelocity with altitude has been computed by Loh [Ref. 7] for 
various values of ballistic coefficients. Figure 16 shows a 
plot of the various deceleration curves. A similar curve can 
be generated from the simulated reentry. The generated curve 
 ЕШІП be the same shape as Figure 16 but since the reentry 
velocity is only about 80 percent of escape velocity, all 
values of dimensionless velocity will be lower. 

In addition, the simulated reentry body undergoes 
ablation and hence changes the value of the ballistic 
coefficient. However, if a nominal ballistic coefficient 


of 1500 is chosen, the value of 





E Jg 50 0016 57) 
NEED 0 8 007-117725 | ( 


can be used, where 9 corresponds to reentry angle. Qualitita- 
tlvely, the simulated reentry curve should lie between the 
0.001 and 0.01 curves and to the left of the curves in 
Figure 16. 

Figure 17 is a graph of the Mach number versus 


altitude for the model. Maximum deceleration occurs at 
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400 


300 


200 


Altitude, 4 ft 


100 





Figure 16. Velocity Variation during Direct Penetration of 
the Earth's Atmosphere. 


(Reproduced from Loh [Ref. 7], p. 110) 
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about 75,000 feet, and the actual impact Mach number is 4.42. 
Based on the above parameters and the similarity in shape 
between curves in Figures 16 and 17, the computer model was 
accepted as a valid representation of the actual reentry 


Brajectory. 


bee TRAJECTORY THROUGH AN ICE CLOUD 
E Model Cloud 

In order to evaluate realistically the effects of 
mS Clouds on the trajectory of a reentry body, a repre- 
sentative ice cloud was needed. 
thickness, particle size and concentration was needed. 
Although information such as this is readily available for 
relatively low altitude rain clouds, useful information 
 Батпітр high altitude ice clouds is very rare. The only 
Buormation found is contained in Table 5, reproduced from 


Heymsfield [Ref. 8], in which aircraft measurements of five 


different cirrus uncinus clouds were made. 


Table 5 


Summary of Cloud Measurements 
(Reproduced from Heymsfield [Ref. 8], p. 800) 


Concen- 
tration 
Aver- 
аме 
| larger 
Sampling than Calculated Calculated 
; Cirrus al- Sampling temperature Ice water 0.15 reflectiv- precipita- 
Flight titude altitude range content* Total mm itv factor* tion rate* 
no. Date Cloud type (km) (km) (°С) (zm m^?) (cm ^3) (cm 73) (mms m7) (mm hr) 
14 *9 July 1970 Ci unc 6.5-8.8 6.5-8.0 —-15 t0 —30 0.39 0.24 — (0.019 17 9 10 0.6 
C(t) 19 Sept 1971 Ci unc 8.5-9.1 8.5-9.1 —36 to —43 0.39 0.24 — 0.025 36 14 1.6 0.7 
1212) 19 Sept 1971 Ci unc 8.9-0.35 8.9-9.35 40 to —44.5 0.25 90.15 — 0.03 14 7 0.97 0.5. 
елер 6 Apr 1972 Ci unc 9.6-11.0 9.6-11.0 —41.7 to —57.6 0.43 0.30 0.40 0.03 29 10 1.0 0.6 
G2) 6 Apr 1972 Ci unc 9.8-10.1 9.8-10.1 —4310 —46 0.30 0.20 0.20 0.04 12 9 0.9 0.6 
етеу б Арг 1972 Ciunc—Cs 9.35-9.8 9,35-9,8 —14:0 —45 0.16 0.075 0.30 0,015 14 4 0.7 0.2 
117 4 Mar 1972 Cs 4.45-8.5 4.45-7.9 —20 to —46 0.067 0.026 — 0.008 4.3 05 70.2 0.1 
E 4 Apr 1972 Cs 7.9-8.9 7.9-8.4 -35 (0 --44.5 0.04 0.025 0.20 0.004 0.6 0.4 0.1 0.025 
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Such information as altitude, 





The table shows that the clouds ranged in altitude 
ШОШ 1 15 то 11 kilometers (14,600 to 36,100 feet) and 
varied in thickness from 0.3 to 3.45 kilometers (980 to 
11,300 feet). Based upon the above figures and the advice 
of Dr. K. L. Davidson, Meteorology Department, Naval 
Postgraduate School, the representative ice cloud was 
selected to be at an altitude of 30,000 feet with a thickness 
of 5,000 feet. Ice particle concentration at that altitude 
nominally was chosen as 0.24 gram per cubic meter which 
corresponds to 4.65 x Ші slug per cubic foot. 

2. Mass Defect 
The mass of the reentry body ablated by collision 


with the ice particles can be expressed by 


Base area x Cloud thickness x Mass Loss Factor 


Mass Defect - SOS 


(58) 


where the base area of the reentry body and the angle of the 
trajectory with respect to the horizontal, A, are computed 

at the proper altitude by the computer program. The Mass 

Loss Factor represents the ratio of the mass ablated to the 
mass of the ice particles which collide with the body. Some 
information for Mass Loss Factor exists for multiple impacts 
without a shock wave ahead of the reentry body, but none could 
be found which included the shock wave. Therefore, several 


nominal values of Mass Loss Factor were selected, and the 





Eu of 10, 20, 350, and 40 was computed. A trajectory 
was calculated for each. 

One assumption inherent in equation (58) is that all 
of the ice particles contained in a tube representing the 
base area of the reentry body passed at the proper angle 
through the cloud thickness do impact on the body. Such an 
assumption was made since this analysis was based on a worst 
case scenario. If only a portion of the ice impacts on the 
body, which is the more likely case, less ablation will occur 
and, therefore, less of an effect will result. 

See Flow Chart 

Using the cloud information contained in Table 5 
and equation (58), the flow chart in Figure 14 was modified 
to include ablation due to ice erosion. Figure 18 describes 


the process. 
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Flow Chart of Trajectory Program with Ice Cloud 





IV. RESULTS, CONCLUSIONS AND RECOMMENDATIONS 


A. RESULTS 


Table 6 compares the difference in range of the various 


trajectories modeled as functions of Mass Loss Factor, which 


is the 
to the 
values 


Figure 


ratio of mass loss of ablation material due to erosion 
mass of ice particles impacting. Also shown are the 
of mass loss due to ice erosion and final Mach number. 


19 shows the same results graphically. 


Table 6 


Comparison .of Trajectories with Ice Erosion 


to Baseline Trajectory 


| 


Mass Loss Factor Final Mach Number Mass Loss Due to Ice (slugs) Difference in Range (feet) | 


0 


0.0947 


0.1895 


0.2842 


0.3789 





B. CONCLUSIONS 


The most interesting result of the analysis is that the 


reentry vehicle with ice erosion traveled farther in all 


cases than did the vehicle without ice erosion. The reason 


is that, despite the significant mass loss due to erosion 


Which tends to decrease the ballistic coefficient, the final 
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nose shape, as seen in Figures 10 and 15, reduces the drag 
coefficient and tends to increase the ballistic coefficient. 
The result is that the two parameters, mass loss and drag 
coefficient, are of opposite sign and tend to cancel. 
“acular Error Probability, CEP, is a constant multi- 


plied by the variance in range, and for a Gaussian distribution 


q = Ig. (597 
Т 


where gi denotes the variance in range caused by the various 
factors that influence changes in range. The assumption can 
be made that another 9i, that is, the variance due to ice 
erosion, can be added to equation (59) and further that the 
variance due to ice erosion is of the order of the change in 
range due to ice erosion. If it also is assumed that the CEP 
15 of the order of the variance, then equation (59) can be 
written 


СЕР = £ CEP,“ = CEP,“ + A Range“ (60) 


1 


where CEP, ‘Sethe original СЕР of the body due to all 
factors excluding ice erosion and A Range is the difference 
in range due to ice erosion. For a C-4 Trident with an 
initial CEP of 1500 feet (from Table 2) and a Mass Loss 
Factor of 40 for a A Range of 15.03 feet, equation (60) 


yields a new CEP of 1500.08 feet. 
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The assumptions in this thesis were made on the basis of 
the best information available in the literature. The analysis 
of the problem based on the assumptions indicates that the 
difference in CEP due to ice erosion is very small indeed. 
Considering ice erosion as the worst case for atmospheric 
water phenomena, the CEP change due to water droplet erosion 
or due to differences in humidity also can be considered to 


be very small. 


C. RECOMMENDATIONS FOR FUTURE WORK 

1. The nose shape after impact with ice particles is 
obviously very important. Experiments which measure the 
effects of multiple collisions of ice on nose shape are 
required. Such experiments must include the effects of a 
bow shock around the nose equivalent to a Mach number of 
about 10. Conceptually, more ice should impact near the 
centerline of the vehicle and less toward the edges, tending 
to make the body more blunt and therefore tending to increase 
tue drag coefficient. 

2. This analysis should be redone using data for a 
carbon-carbon nose and new data for ice erosion found in 
part 1 above. The mass loss for a carbon-carbon nose tip 
would tend to be far less than for the silaceous nose tip 
which is modeled in this thesis. Hence, mass loss due to 
general erosion would tend to have less effect than in this 


analysis. 
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It may be, however, that the mass loss due to ice erosion 
for the carbon-carbon nose is of the same order of magnitude 
as the mass loss due to ice erosion for the silaceous nose. 
If so, the effects of ice erosion may be more important in 
the carbon-carbon nose than in the silaceous nose because 
the effects of mass loss due to ice erosion were masked in 
the silaceous nose by the substantial mass loss occurring 


throughout the reentry. 
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APPENDIX À 


POLYNOMIAL REPRESENTATION OF NOSE SHAPES 


The following pages contain graphs of the various nose 
shapes. Each graph contains a plot of the data points taken 
from the respective nose shape curves in Figure 10 with the 
aid of the HP-9874A digitizer and a plot of the polynomials 
which describe the data points. The polynomials also are 
displayed on the graph; the fifth-order polynomial applies 
ШЕ? < X < 2.4 and the first-order polynomial applies for 
E 2.4. 

In the analysis, the polynomial was used to describe the 
first 2.4 inches of the nose shape. The linear approximation 
Was used to describe the nose shape from 2.4 inches aft of 
the nose tip to about 20 inches aft of the nose tip. The 
linear approximation to Figure A-1 was used to describe the 
nose shape from about 20 inches aft of the nose tip (where 
the two linear portions intersect) to the end of the reentry 
body for all nose shapes. 

The profiles for nose shapes with ice erosion use the 
same equations to describe the shape. The coefficients are 
modified by substitution of a modified altitude into the 


equations for the coefficients displayed in Appendix B. 


Z 1 
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DEDENDIX B 


ГОК 22 REPRESENTATION OF COEFFICIENTS 


The following pages contain graphs of the various nose 
shape coefficients defined in equations (49) through (52). 
The first six graphs correspond to the coefficients used in 
Stations (49) and (51) and apply for 0 < x < 2.4 inches 
where x is centerline distance from the nose tip. The last 
two graphs correspond to coefficients used in equations (50) 


ШІ (52) апа агріу for x > 2.4 inches. 
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APPENDIX C 


POLYNOMIAL REPRESENTATION OF NOSE TIP RECESSION 


The following graph shows the evaluation of nose tip 
recession, the distance the nose tip ablates from the 
original location, as shown in Figure 10. As for the 
coefficients described in Appendix B, the nose tip recession 
was divided into two portions. А sixth-order polynomial 
MeSetibes the nose tip recession below 90,000 feet апа а 


linear model describes the recession above 90,000 feet. 
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APPENDIX D 


COMPUTER PROGRAM EXAMPLE 


Appendix D is a sample of the computer program used to 
conduct the trajectory analysis. The computations were 
conducted on an HP-9845A computer which uses BASIC computer 
language. In general, the program is self-explanatory in 
that the names of variables are used whenever possible. 

Line numbers beginning with 1XXX describe the initial 
conditions. The counter K is zero until zero altitude 15 
reached and then is equal to one. The counter L is zero 
during the first iteration and allows the establishment of 
initial volume, mass, and base area of the reentry body. 
After the first iteration of nose shape, L is equal to 1. 
The counter Z is zero until the altitude of the cloud is 
reached. The value of Z is changed to one during the 
lteration to find the mass loss due to the impact of the 
ice in the cloud. After the iteration is complete, the 
value of Z is two. If Z is set equal to two initially, the 
baseline trajectory with no ice ablation will be calculated. 

Line numbers beginning with 2XXX calculate the new 
position of the body as it descends. Line numbers beginning 
With 4XXX calculate all of the major parameters needed for 
the calculation of pressure coefficients and the nose shape 
coefficients. Line numbers beginning with 5XXX calculate 


the actual nose shape, pressure coefficient, drag, volume, 
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etc. Line numbers beginning with 6XXX calculate the initial 
volume, mass and base area and values of drag coefficient 
and ballistic coefficient during the entire reentry. Line 
numbers beginning with 7XXX print the final values obtained 
during reentry and stop the program. Line numbers beginning 
with 8XXX calculate the mass ablated during passage through 
the ice cloud. Line numbers beginning with 9XXX calculate 
the final range and set the final altitude to zero so that 


calculation of other parameters at zero altitude is possible. 
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